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The potential energy surface (PES) of dissociative adsorp- 
tion of H2 on Pd(lOO) is investigated using density functional 
theory and the full-potential linear augmented plane wave 
(FP-LAPW) method. Several dissociation pathways are iden- 
tified which have a vanishing energy barrier. A pronounced 
dependence of the potential energy on "cartwheel" rotations 
of the molecular axis is found. The calculated PES shows no 
indication of the presence of a precursor state in front of the 
surface. Both results indicate that steering effects determine 
the observed decrease of the sticking coefficient at low ener- 
gies of the H2 molecules. We show that the topology of the 
PES is related to the dependence of the covalent H(s)-Pd(d) 
interactions on the orientation of the H2 molecule. 

PACS: 68.10Jy,68.35Dv,73.20At,82.65Jv 



I. INTRODUCTION 

Hydrogen dissociation on metal surfaces has become 
one of the benchmark systems for the understanding of 
dissociative adsorption processes. The H2 / Pd system 
is particularly interesting as an example of non-activated 
dissociation over a transition metal surface. The disso- 
ciation process has been studied using molecular beam 
adsorption experimentsu as weW as state resolved time=cj£- 
flight measurements of desorbing hydrogen moleculesa~El. 

Hydrogen molecules dissociate spontaneously at 
Pd(lOO), i.e. molecules with low kinetic energy disso- 
ciatively adsorb at|-.Ed(100) surfaces with a large initial 
sticking coefficienttlQ. The preference of non-activated 
dissociation pathways for low kinetic energies of the H2 
molecules is also consistent with the independence of the 
initial sticking coefficient on incident angle 8 as well as 
with the measured Maxwell-Boltzmann velocity distribu- 
tion of desorbing H2 molecules which correspond^ to the 
surface temperature of the palladium substrata3fl. 

There are, however, clear indications that activated 
pathways exist as well. After a decrease for beam ener- 
gies below K, 0.15 e.V the initial sticking coefficient con- 
tinuously increasescl. At those larger kinetic energies the 
dependence of the initial sticking coefhcient on the in- 
cident angle of the molecules changes into a cos"~^(8) 
behavior with n increasing up to n « 1.6 for beam en- 
ergies of 0.4 eV. State resolved desorption experiments 
show that the occupancy of the first vibrational excita- 
tion is significantly higher than that corresponding to a 
gas of hydrogen jintilecules in equilibrium with the sur- 
face temperaturaa'DH. This finding is called vibrational 



heating. On the other hand, the rotational degrees of 
freedom are occupied less than what the surface temper- 
ature would imply (rotational cooling)Q. In order to ex- 
plain the vibrational heating, the values for n > 1 in the 
angular distribution, as well as the increase of the stick- 
ing coefficient with the kinetic energy of the H2 molecules 
it has been found necessary in dynamic studies on two- 
dimensional model potential-energy surfaces (PES) to as- 
sume the presence of a small (about 0.1 eV) barrier in the 
exit channel^ i.e. in the region of the reaction path™ 
where the H-H bond is already significantly strctchcdH 
Such an energy barrier, however, yields an additional ki- 
netic energy of thermally desorbing molecules which has 
not been measured experimentallyo. 

There has been a long debate about the origin of the 
initial decrease of the initial sticking coefficient with in- 
creasing beam energies. Such a decrease and the observed 
dependence of the initial sticking coefficient on the angle 
of incidence with a coefficient n < \ have been usually 
assigned to the presence of a weak, metastable molec- 
ular bound state of the H2 moloGultj-jB, large distance 
from the surface (precursor statc)lil-t2HiirO. The presence 
of a precursor state has also been assumed in order to 
explain the slow decrease of the sticking coefficient with 
hydrogen coverages. An alternative explanation relates 
the initial decrease to the presence of both acti^pied. and 
non-activated dissociation pathways on the PESEjlj. At 
low kinetic energies molecules with unfavorable orien- 
tations may reorientate during the dissociation process 
and follow a non-activated pathway but with increasing 
kinetic energy this steering is more and more hindered 
and the number of molecules proceeding along activated 
pathways increases. 

The fundamental understanding of the dissociative ad- 
sorption of hydrogen rests on the knowledge of the PES 
which, as will be stressed below, depends rather crui- 
cially on all six coordinates of the two hydrogen atoms 
of the H2 molecule. Experimental data have been inter- 
preted in terms of simple, often two-dimensional model 
PES as discussed above. Very recently^ccurate den- 
sity functional calculation (DFT) of PESBTlapd Mgh- 
dimensional quantum-dynamical simulationaU'Ej Ej of 
H2 dissociation have become available thus opening the 
way for complete ab inito investigations of dissociative 
adsorption and associative desorption processes. Previ- 
ous DFT studies of the H2 interaction with metal sur- 
faces concentrated mainly on system where the dissoci- 
ation is activated. There are a few studies of the inter- 
action of H2 molecules in front of surfaces of transition 
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metals with a partly filled d-statesE-jOta. In the recent 
work of White and Birdll3 investigating hydrogen disso- 
ciation on W(IOO) the presence of non-activated dissoci- 
ation pathways was confirmed but in dependence of the 
center-of-mass position of the molecule in the surface unit 
cell even for the most favorable orientation of the molec- 
ular axis activated dissociation has been obtained as well. 
Molecules with its axis perpendicular to the surface have 
been found to experience no attractive interaction. Sum- 
marizing their results, White and Bird expect molecular 
steering to play an important role in the dissociation dy- 
namics at low incident energies. 

In this paper we present and discuss detailed calcula- 
tions of the potential energy surface for hydrogen disso- 
ciation on Pd(lOO). Here we analyze the physical prop- 
erties behind the topology of the PES. We note that the 
results are now sufficiently complete for a six-dimensional 
quantum-dynamical simulation of the dissociation of II2 
molecules on Pd(100)E3. 

II. CALCULATION METHOD 

The calculations have been performed using density- 
functional theory, within the generalized gradient approx- 
imation (GGA)cJ and the full potential linear augmented 
plane wave method (FP-LAPW)cSEa. The application 
of the GGA is of crucial importa nce , fo r the calculation 
of PES of hydrogen dissociationi3'E3EZl,but some impor- 
tant questions concerning the, accuracy of PES calculated 
applying the GGA remainES. The good agreement of 
experimental data with the results of high-dimensional 
quantum-dynamical simulation of II2 dissociation based 
on the calculated PES ior H2/Cu(111)EI and the PES 
presented in this papeiE3 shows that the GGA repro- 
duces the important features of these PES with sufficient 
accuracy. 

The FP-LAPW wave functions in the interstitial re- 
gion are represented using a plane wave expansion up to 
^'cut=ll Ry, and due to the small radius of the muffin-tin 
sphere around H atoms (r^^j^ — 0.37 A) it is necessary to 
take into account plane jssiaves up to £^cut = 169 Ry for the 
potential representationcJ. Inside the muffin-tin spheres 
the wave functions are expanded in spherical harmon- 
ics with Zmax = 10, and non-spherical components of the 
density and potential are included up to Ima.^ = 3. For 
the k-integration we used 64 uniformly spaced points in 
the two-dimensional Brillouin zone corresponding to the 
c(2x2) surface unit cell. All calculations were performed 
non-relativistically. 

In the calculations we used a supercell geometry. The 
metal substrate is modeled by three layers slabs which 
are separated by a 10 A thick vacuum region. Hydrogen 
molecules are placed at both sides of the slab. We used 
a rigid substrate because due to the large mass ratio be- 
tween Pd and H the substrate will not change during the 
scattering event. The slab thickness of only three layers 



has been chosen with respect to the large computational 
demand of the evaluation of the PES. The application 
of such small layer thickness for the Pd(lOO) substrate 
gives an acceptable accuracy of the PES because due to 
the large density of states at the Fermi energy perturba- 
tions in the electronic structure induced by the impinging 
H2 molecule are screened very effectively. We confirmed 
this by increasing the layer thickness of the substrate to 
5 layers and comparing the total energy with respect to 
the separated H2 and Pd(lOO) surface at various points 
of the reaction pathways. We found that in the region of 
the PES where the H2 molecule approaches the surface 
and starts to dissociate the differences of the potential 
energy are less than 0.05 eV per H2 molecule. In partic- 
ular, the character of the dissociation pathway, whether 
it is activated or non-activated, is not influenced. This 
part of the PES determipes the dynamics dissociation of 
the hydrogen molcculcE3. The dependence of the adsorp- 
tion energy on the layer thickness is more pronounced an 
the region of the PES corresponding to adsorbed H atoms 
and strong H-Pd bonds. Here, the energy changes with 
increasing layer thickness up to 0.12 eV per H2 molecule. 

The application of such a small layer thickness for the 
Pd(lOO) substrate is possible because due to the large 
density of states at the Fermi energy perturbations in the 
surface electronic structure induced by the impinging H2 
molecule are screened very effectively. We confirmed this 
by increasing the layer thickness of the substrate to 5 lay- 
ers and comparing the total energy with respect to the 
separated and Pd(lOO) surface to that of the three layer 
calculations at various points of the reaction pathways. 
We found that the differences in the potential energies 
are below 0.12 eV and that the changes are considerably 
smaller for larger heights of the molecule above the sur- 
face. For the PES calculation we employed a supercell 
with a c(2x2) surface structure. However, in the case of 
the geometry where the H2 hits an on-top site and disso- 
ciates into two hollow positions H atoms belonging to H2 
molecules in different surface unit cells within a c(2x2) 
structure would come too close during the dissociation 
process and the unit cell has been enlarged to a p(2x2) 
structure. The convergence of our results was tested us- 
ing a p(2x2) surface unit cell, a plane- wave cut-off for 
the wavefunctions expansion of £'cut=13 Ry and a larger 
H muffin tin radius r^r^ = 0.48 A in the exit channel re- 
gion. The resulting changes of the potential energy were 
again found to be less than 0.08 eV per H2 molecule. 

The potential energy of a hydrogen molecule as it is 
used throughout the paper is defined as the DFT-GGA 
total energy. The energy zero is the energy of the ge- 
ometry difference where the molecule is sufficiently far 
away from the surface, {Z= 3.7 A), such that there is 
practically no interaction between the molecule and the 
surface. Zero point corrections and other vibrations are 
not included in the PES. 



2 



III. RESULTS 

First, we briefly summarize the properties of bulk 
Pd, tlaG|-eJean Pd(lOO) suface, and adsorbed hydrogen 
atomgH3c3. The calculated lattice constant of Pd calcu- 
lated by the DFT-GGA is 4.03 A. This value is 3 % larger 
than the lattice constan obtained within |tlje LDA, and 
3.5 % larger than the experimental valuecS. The clean 
Pd(lOO) surface shows a slight inward relaxation of the 
topmost Pd layer of —0.9 % of the bulk interlayer dis- 
tance. This value is close to the r|eIaxation of —0.6 % 
obtained by DFT-LDA calculationcJ. The equilibrium 
position of the hydrogen is the surface hollow site with 
a small adsorption height Zo=0.1 AEj. At the bridge 
site the adsorption energy is about 0.3 eV less favorable 
and the adsorption height increases to Zo= 1.3 A. Ad- 
sorption at the on-toQjDOsition is even unstable against 
associative desorptiorO. Obviously, this pronounced de- 
pendence of the H - Pd bond strength on the geometry 
of the H adatom effects the potential energy during the 
dissociation process as will be discussed below. 

Assuming that the positions of the substrate metal 
atoms can be kept fixed the PES of hydrogen dissociation 
depends on the six coordinates defining the position of 
the two hydrogen atoms of the H2 molecule. Figure ^ 
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FIG. 1. Coordinate system used to describe the orientation 
of the H2 molecule with respect to the surface Pd atoms. The 
height of the H-H center of mass above the surface is Z, the 
H-H distance is dn-H, and the angle of the molecular axis with 
the surface normal is 9. 



licopter rotation). The distance between the hydrogen 
atoms is denoted by dn-H- 

In order to investigate possible non-activated pathways 
Fig. I gives 
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FIG. 2. Cut through the six-dimensional potential energy 
surface (PES) of a H2 molecule in front of Pd(lOO). We dis- 
play "elbow plots" where Z is the height of the H2 center 
of mass over the surface, and dn-H is the distance between 
the two hydrogen atoms. The cut is defined by keeping the 
molecule parallel to the surface at an azimuthal orientation 
shown in the inset. The units of the potential energy are eV 
and the interval between adjacent contour lines is 0.1 eV. 

the PES as functions of H-H distance, dn-H, and the 
molecular center-of-mass distance from the surface, Z, 
for fixed, characteristic positions of the molecules center- 
of-mass in the surface unit cell: the bridge {X=0, Y=0.5 
dpd-Pd, (/3=0°))(a), the hollow (b) {X=Y=0.5 dpd-Pd, 
(p=0°) and the on-top {X^Y^O, (^=45°) (c) site. The 
molecular axis is kept parallel to the surface {6=90°). 
For all geometries investigated in Fig. || this 9=90° ori- 
entation is the energetically most favorable one (see also 
Fig. 1^ below) . The dissociation over the 



shows the coordinate system used to describe the po- 
sition and orientation of the H2 molecule. The center of 
mass position is given by the three cartesian coordinates 
{X, Y, Z) for which the origin is at the center of a sur- 
face Pd atom. The two rotational degrees of freedom are 
given by the angle of the molecular axis with the Z-axis 
(cartwheel rotation) and with the X-direction Lp (he- 
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FIG. 3. Dependence of the potential energy of an H2 
molecule on the on the angle Q of its molecular axis with 
the surface normal. The geometry for 61=90° corresponds to 
Fig. 2a. The molecule is in the entrance channel at (Z=1.6 A, 
ciH-H=0.8 A) (open circles) or in the exit channel at {Z—1.0 
A, dH-H=l.l A) (filled circles) 

bridge (Fig ^) and over the hollow site (Fig. ||b) is 
non-activated and proceeds with a continuous gain of en- 
ergy without a hampering energy barrier. The dissoci- 
ation of an hydrogen molecule over the on-top site with 
fixed coordinates (^=0, Y=0) is activated. The barrier 
has an energy of about 0.15 eV and is situated in the 
exit channel region. Figure ^ shows that non-activated 
dissociation of H2 molecules on Pd(lOO) is possible over 
several non-activated pathways including a large region 
of center-of-mass positions of the hydrogen molecule in 
the surface unit cell. 

The two-dimensional cuts through the PES in Fig. |^ 
refer to the situation that the orientation of the molecular 
axis is fixed. In reality the impinging hydrogen molecule 
rotates. In the following we discuss the dependence of 
the potential energy on the orientation of the molecular 
axis, i.e. on the angles 9 and (p. 

The change of the orientations the molecular axis with 
respect to the surface normal, 9, and away from 6—90°, 
implies a significant increase of the potential energy. Fig- 
ure m displays the change of the potential energy of a 
H2 molecule in dependence of the angle 9 of the molec- 
ular axis with the surface normal for two positions in 
the PES of Fig. all other parameters beeing fixed. 
For both distances, Z=1.6 A and Z= 1 A the potential 
energy increases like cos^(0). Fig. || shows that with de- 
creasing distance to the surface a minimum in the depen- 
dence of the potential energy with angle 9 with increas- 
ing depth develops. The "cartwheel" rotations are more 
and more hindered. There is, however, a large difference 
between the "entrance" channel region, where the hydro- 
gen molecule approaches the surface and the II2 bond is 
still intact and the "exit" channel region where the H-H 
bonds is stretched and the dissociation of the molecule 



starts. In the entrance channel, at larger heights Z (see 
open circles in Fig. |^), there is no or only a small repul- 
sive interaction of the molecule with the surface even for 
a configuration with the molecular axis perpendicular to 
the surface 9—0°. Thus, the hydrogen molecules may ap- 
proach the Pd(lOO) surface without an appreciable loss in 
energy and independent of their initial rotational state. 
In the exit channel where bonds of the hydrogen atoms 
with the surface form the potential energy depends sensi- 
tively on the orientation 9 (see filled circles in Fig. |3|) and 
the interaction of the molecule with the surface becomes 
repulsive already for small deviations from 0=90°. These 
conclusion are confirmed for the extreme case of a hydro- 
gen molecule hitting the surface at the hollow site with 
its axis fixed perpendicular to the surface. We find that 
in the entrance channel the downmost hydrogen atom 
may approach the surface plane without a substantially 
barrier but for fixed orientation of the molecular axis at 
9=0° the potential energy increases drastically if the H- 
H bond is stretched and a nearly classical hard wall po- 
tential blocks the exit channel. Without any steering of 
the molecular axis during the dissociation non-activated 
dissociation pathways are possible only for molecules en- 
tering the exit channel region within a narrow region of 
orientations 9 of the molecular axis and low initial stick- 
ing coefficient would result. 

We expect, however, that the steering of the molecular 
axis toward energetically favorable orientations is rather 
effective. At low kinetic energies the rotational period 
of the molecule is comparable to the time the molecule 
spends at distances Z corresponding to the end of the 
entrance channel. A molecule with a kinetic energy of 
£;kin=0.05 eV travels 0.6 A per 0.03 ps and when it is in 
the j=2 rotational state it rotates during this time (or 
pathway) by about 90°. The large increase of the po- 
tential energy with angle 9 in Fig. || at the end of the 
entrance channel corresponds to large forces aligning the 
molecular axis towards orientations where the dissocia- 
tion becomes energetically favorable ("steering effect"). 
With iacreasing beam energy in the adsorption exper- 
imentsEl both the kinetic and rotational energy of the 
molecules increases. This shortens the time the molecule 
spends at the end of the entrance channel and it becomes 
more likely that the molecule is reflected before the axis 
is aligned. An increase of the rotational energy implies 
that states with shorter rotational periods are occupied 
and those molecules turn out of the favorable orienta- 
tion again more easely. As a consequence of th&ireduced 
steering the inital sticking coefficient reducest[c3. 

In the high symmetry configurations shown in Fig. ^ 
the most favorable orientation of the II2 molecule is 
that where the molecular axis is parallel to the surface 
(6'=90°). For other center-of-mass positions and orien- 
tations (p of the molecule the minimum of the potential 
energy with respect to the angle 9 deviates from 6'=90°. 
As an example, if the center of mass of the molecule is 
between the hollow and bridge position (X=0.25 dpd.pd, 
y=0.5 rfpd-Pd, V=0°) the optimal angle at a height 
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Z=l.l A is 9=65°, i.e. the H atom pointing toward 
to hollow site is closer to the surface than that point- 
ing toward the bridge site. This tilting of the molecular 
axis is consistent with the different adsorption height of 
hydrogen adatoms at the hollow and bridge adsorption 
site. The potential energy is -0.3 eV, a value between 
that found along the reaction pathway for a center of 
mass-position over the bridge (Fig. ||a) and hollow site 
(Fig. I^b) at this height. Generally, the energetically most 
favorable orientation of the molecular axis, 0, optimizes 
the H-surface bond strength of both H atoms of the H2 
molecule. 

In the case of "helicopter" rotations a pronounced de- 
creases of the potential energy is found for geometries 
where H and Pd atoms get closer than a typical H-Pd 
bond length (1.7 A). These geometrical constraints af- 
fect different positions in the surface unit cell differently. 
Dissociation over the bridge site is non-activated in the 
geometry shown in Fig. ga (i^=0°). In the case that 
the hydrogen atoms are oriented towards the Pd atoms 
{ifi=90°) helicopter rotations are not hindered only at 
large distances Z. For this configuration a steep increase 
of the potential energy is found already at the end of 
the entrance channel (Z=1.2A, dH-H=0.9 A). The H-Pd 
distance is 1.54 A and the potential energy for H atoms 
orientated towards the Pd atoms is i?pot=+0.25 eV, i.e. 
0.57 eV higher than in geometry with the H atoms ori- 
ented toward the hollow sites (see Fig ^). In the case 
that the center of mass position of the impinging molecule 
is over the hollow site helicopter rotations are not hin- 
dered. A substantial increase of the potential energy for 
the orientation of the molecule with the H atoms oriented 
toward the surface Pd atoms ((^=45°) is found only in the 
exit channel where the H-H bond is already broken and 
both H atoms may easily reorientate separately. In dif- 
ference to the "cartwheel" rotations the importance of 
the steering of the impinging hydrogen molecule into fa- 
vorable orientations ip of the molecular axis will depend 
on the position of the molecule in the surface unit cell. 

Metastable, weakly bound molecular states of H2 
which may act as precursor states to the dissociation 
of the molecule are commonly assumed to arise due to 
an interplay the repulsive Pauli repulsion between the 
closed shell electronic configuration of H2 and the tails 
of the metal surface electronic states and the attractive 
van der Waals interactioncl. In our calculations there 
is no indication of local minimum of the potential en- 
ergy related to such types of precursor states at large 
distances in front of the Pd(lOO) surface. Although the 
long-range tail of the van der Waals interaction related 
to polarization fluctuations are not described within the 
GGA recent calculations-ipf the PES of He atoms in front 
of Rh(lOO) have shownEj, that the GGA is cabable to 
adequately reproduce weakly bound states of adsorbates 
with a closed shell electronic configuration. The absence 
of a considerable repulsive interaction of the impinging 
H2 molecule with the surface at larger distances Z even 
for unfavorable orientations of the molecular axis and the 



large decrease of the potential energy already at distances 
Z «2 A found in Fig. || make this mechanism of the for- 
mation of molecular bound states in front of the metal 
surface rather improbable. 

Local minima in a restricted subspace of the high- 
dimensional PES, may arise not only at large distances 
in front of the surface due to the interplay of the Pauli 
repulsion and attractive van-der Waals interaction but 
also at distances Z where the direct H-metal interaction 
dominates. They are a consequence of the pronounced 
dependence of the H-metal bond on the orientation of 
the molecule and the H-surface-metal-atom distances. 
For an example, the dissociation pathway on top of a 
Pd surface atom (Fig. ^) shows a minimum of -Epot m 
the entrance channel for fixed center-of-mass position of 
the molecule. 

A more detailed investigation of the dissociation over 
the on-top site shows that the minimum of potential en- 
ergy seen in Fig. Ic is not a local minimum in the PES. 
The molecule can follow a purely attractive pathway if 
its center-of-mass position is able to relax in order to op- 
timize the H-Pd bonding. Figure |4| shows the variation 
of the potential energy for 
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FIG. 4. Dependence of the potential energy of H2 on 
its center-of-mass-position along a path R={X=~t, Y=t, 
^=constant) for three heights of the molecule. The molecu- 
lar axis is parallel to the surface (^=90°) and the angle is 
fixed at <^=45°. The point t=Q corresponds to the geometry 
in Fig. 2c. 

different heights Z in the case that the center- 
of-mass coordinate is changed along a path, i, 
i?=(X=— i,F=i,Z=constaiit). At large distances Z > 
1.8 A a position on-top of the surface Pd atom {t=0) is 
most favorable. At smaller heights the optimal reaction 
pathway involves a motion of the center-of-mass position 
towards the hollow sites. This geometry for lower Z is 
similar to the energetically most favorable ocientation of 
a H2 molecule in front of a Rh(lOO) surfaceO. 

It is interesting to-compare the PES of hydrogen dis- 
sociation on W(100)lla with that on Pd(lOO). Both sub- 
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strates are transition metals but tungsten and palla- 
dium crystalize into different lattice types (W is bee, 
and Pd is fee). Furthermore, W is in the middle of 
the 5d and Pd at the end of the Ad transition metal 
series. The PpS of hydrogen dissociation on Pd(lOO) 
and on W(100)t3, however, is qualitatively similar. Both 
PES exhibit several non-activated dissociation pathways. 
The pronounced dependence of the potential energy on 
the orientation of the molecule with repect to the sur- 
face metal atoms, especially on "cartwheel" rotations, is 
found also on W(IOO). There are also characteristic dif- 
ferences between both PES. This is already documented 
in the different equilibrium adsorption site, which is at 
the bridge position on W(IOO) in contrast to the hol- 
low site on Pd(lOO). The geometry of the non-activated 
dissociation pathways on both surfaces is different. For 
example the dissociation with fixed center of mass posi- 
tion over the hollow site has an energy barrier of 0.3 eV 
on W(IOO) but no barrier on Pd(lOO). These differences 
will influence details of the dynamics of the hydrogen 
molecules on the PES. The calculations confirm, how- 
ever, that the central features of the PES of hydrogen 
dissociation, e.g. the presence of non-activated dissocia- 
tion pathways, are determined by the presence of a partly 
filled d-band in the electronic structure of both transition 
metal substrates. 



IV. BREAKING OF THE MOLECULAR BOND 
AND MAKING OF NEW ADATOM-SURFACE 
BONDS 

The large reactivity of transition metals for hydrogen 
dissociation given by the presence of dissociation path- 
ways with no or only sipj Jl . en oMy barriers is actuated 
by a partly filled d-bandEa't§E3~E3. The mechanism how 
the d-band may affect the energetics of hydrogen disso- 
ciation has been interpreted— in different ways. Within 
the Harris-Andersson mode£3 the repulsive interaction 
and the formation of barriers for H2 dissociation on sur- 
faces of simple and noble metals is, based on the work 
of Zaremba and KohnEZi, related the Pauli repulsion be- 
tween the closed shell configuration of the H2 molecule 
and the tails of the s,p- waves of the surface electronic 
structure. At transition metal surfaces the Pauli repul- 
sion is suppressed by a low-energy occupation redistribu- 
tion of electrons out of the s,p-tails into the large density 
of holes in the d-band just above the Fermi energy. Thus 
there is no direct interaction with the metal d-states and 
the dissociation barriers are governed by the interaction 
at large distances to the surface. The analysis of PES 
of hydrogen dissociation on various metal surfaces has 
shown that the top of the energy barriers occurs in the 
exit channel, i.e. at a point were the H-H interaction is 
significantly weakened_and.and strong H-metal bonds are 
already beeing formedll3'0. Those results, which are not 
consistent with the Harris-Andersson modelE3, hav^ibewi 
consistently explained within a reactivity modelllafiJ^ES 



related to the frontier orbital conceptS'S. Superimposed 
to the interaction of the hydrogen molecule with the s,p- 
states of the metal, which gives rise to energy barriers 
found at simple metals, there is a contribution due to 
the direct interaction of H2 with the metal d-states. It is 
repulsive as long as both the bonding and the antibond- 
ing states of the interaction of the H2 CTg-orbital with 
the metal d-bands are occupied but becomes attractive 
in the case that the antibonding states get depleted. For 
transition metals with a high density of d-states just at 
and below the Fermi level the states which are antibond- 
ing with respect to the adsorbate surface interaction are 
pushed above the Fermi level. Thus these antibonding 
states remain empty level already in the case of a weak 
H2-surface interaction and a small or absent dissociation 
barrier results, 

Feibelmantj calculated the dependence of the poten- 
tial energy of a hydrogen molecule in front of a Rh(lOO) 
surface in dependence of the orientation of the molecule 
with respect to the surface metal atoms for fixed heights, 
Z. He found, that the potential energy variations are gov- 
erned by the orientation dependence of the bond strength 
of the H atoms and the metal d-state. The most favorable 
orientation of the H2 molecule is not that of the high- 
est symmetry but that optimizing the H-surface bond 
strength. 

Our results confirm the conclusion that the topology 
of the PES is determined by the direct interaction of the 
hydrogen molecule with the electronic states of the metal 
surface. The PES refiects the subtle balance between the 
weakening of the H-H bond in the molecule and the en- 
ergy gain connected wit h th e rehybridizing into bonds 
with the metal substrateli3E3. Therefore, it is not un- 
expected that a correlation exists between the energy at 
a given height, Z, along the reaction path and the en- 
ergy and optimal adsorption height of a hydrogen adsor- 
bate at the position of the atoms of the H2 molecule. 
The end of the entrance channel regions in Fig. I and, 
hence, the height Z of a hydrogen molecule at which a 
significant stretching of the intermolecular H-H distance 
occurs, correlates with the adsorption height of a sin- 
gleJaydrogen adatom at the center of mass position (X, 
y)E3. At large distances, Z > 1.9 A, center-of-mass po- 
sitions of the hydrogen molecule on-top of a Pd atom are 
energetically most favorable and center-of-mass position 
over the hollow site have the lowest gain in potential en- 
ergy. The direct dissociation with the center-of-mass of 
the molecule over the bridge site into two neighboring 
hollow sites (dH-H= 2.85 A) in Fig. ^ is connected with 
a large, gain in potential energy and a small curvature in 
the "elbow" plot which is consistent with the large differ- 
ence in adsorptian height and energy between bridge and 
hollow positionEj. In the case that the molecule impinges 
at the hollow site (Fig. ^d) the situation is reversed. The 
molecule first closely approaches the surface but during 
the dissociation process the hydrogen atoms move up- 
ward to the energetically less favorable bridge positions. 
Consequently, although the dissociation pathway shown 
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in Fig. lb is non-activated as well, the energy gain during 
dissociation process of the hydrogen molecule is smaller 
and the curvature in the "elbow" plot is larger. The 
impact position of the molecule over the on-top site of 
a surface Pd atom (Fig |^) is energetically most favor- 
able for larger Z. However, the H-Pd interaction on-top 
of a metal atom is not strong enough to break the H-H 
bond without the formation of an energy barrier which is 
in correspondence with the calculated instability of hy- 
drogea adatoms against associative desorption at on-top 
sitesc3. 

A more detailed picture of the interaction of a hydro- 
gen molecule with the metal surface may be obtained 
analyzing the electron density induced by a hydrogen 
molecule in front of the surface. In Fig. || we show the 



(a) (b) 




FIG. 5. Difference of the charge density between a hydro- 
gen molecule in front of Pd(lOO) at a height Z=1.8 A and that 
of the superposition of the charge density of the clean surface 
and of free hydrogen molecules. The center-of-mass position 
of the molecule is over the on-top site (a) and the hollow site 
(b). The cut plane is perpenticular to the surface and con- 
tains the II2 molecule. The units for the charge density are 
10~^ A"'' and A for the height above the surface. 

difference of the electron density between a hydrogen 
molecule in front of Pd(lOO) at a height Z=1.8 A and 
that of the superposition of the density of the clean sur- 
face and of free hydrogen molecules. Two different con- 
figurations are considered. In Fig. ^ja the H2 molecule is 
over the on-top site and the potential energy has already 
a large negative value of -0.28 eV and corresponds to the 
minimum in Fig. In comparison Fig ^jb shows the sit- 
uation where the H2 molecule is over the hollow site and 
the potential energy is close to zero. Figure |^a clearly 
visualizes that the Il2-metal surface interaction is con- 
nected with a large polarization of both the H2 molecule 
and the metal valence electrons. The drop in the po- 
tential energy in Fig. ||c is caused by a direct interaction 
between the dg state of the II2 molecule and the — r^- 
orbital of the surface Pd atom whickli|is a large density 
of states just below the Fermi levelE3C3. There is an oc- 
cupation redistribution inside the d-shell which depletes 
the 3z^ — r^-orbital thus reducing the Pauli repulsion. 
Furthermore, holes are introduced into the dg state weak- 
ening the H-H bond. In addition an increase of electron 



density around the center of the H atoms and a more de- 
localised accumulation of charge density around the Pd 
atom occurs corresponding to an interaction of the anti- 
bonding cr*-orbital of the hydrogen molecule with metal 
states just above the Fermi level. In the case that the H2 
molecule is shifted from the on-top towards the hollow 
site in Fig. ^ the interaction with surface Pd d-states 
having a large density of states just below the Fermi en- 
ergy is much weaker. The charge density induced around 
the hydrogen molecule is now more similar to the orthog- 
onalization hole found fop-hydrogen dissociation at the 
s,p-metal surface Al(110)c3. But also in this case there 
is a considerable depletion of the 3z^ — r^-orbitals of the 
surface Pd atoms and an internal occupation conversion 
inside the d-shell. Due to the weaker polarization of the 
d-orbitals the potential energy is increased in comparison 
to the position over the on-top site but also in this case 
no barrier forms. It is interesting to note, that the po- 
larization pattern in Fig. |^ is close to that deduced from 
the electronic properties of the non-interacting Pd(lQI]L) 
surface using the reactivity function introduced in Ref.Lj. 
The induced charge density shown in Fig. ^ is consistent 
with the reactivity model developed in Ref. |l8|,^ 

V. CONCLUSION 

We have calculated the PES of hydrogen dissociation 
on the clean Pd(lOO) surface using density functional the- 
ory and treating exchange-correlation in the Generalized 
Gradient Approximation (GGA)c3. The results confirm 
the presence of several, non-activated dissociation path- 
ways. It is shown that "cartwheel" rotations of the molec- 
ular axis out of the favorable orientation parallel to the 
surface are connected with a large increase of the po- 
tential energy, thus confining non-activated dissociation 
pathways to a small angular region of molecular axis ori- 
entations nearly parallel to the surface. "Helicopter" ro- 
tations, on the other hand, are hindered mainly because 
of steric restrictions if the H-Pd distance becomes shorter 
than a typical H-Pd bonding length. The calculations 
gave no evidence of a weak molecular bound state, or 
precursor state, of the H2 molecule at large distances in 
front of the surface. There is nearly no repulsion of H2 
molecules in the entrance channel of the reaction path- 
ways even in the case of the activated adsorption path- 
ways and thus no appreciable energy is required for the 
hydrogen molecule to come close enough to the surface 
in order to form hydrogen-metal bonds. The observed 
pronounced differences in the potential energy in depen- 
dence of the geometry of the molecule with respect to the 
surface is governed by the orientation and bond length 
dependence of the covalent interaction of the hydrogen 
atoms mainly with the different Pd d-bands. Summa- 
rizing these results it is concluded that the dynamics of 
the hydrogen molecule on the high-dimensional PES can- 
not be understood using low-dimensional model surfaces 
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but require a quantum dynamic simulation including all 
degrees of freedom of the hydrogen molecule. Those cal- 
culations have been performed usirg-the calculated PES 
presented in this paper as an inputE3. 
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